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Abstract: Exploiting i?-meson decays for Standard Model tests and beyond requires 
a precise understanding of the strong final-state interactions that can be provided model- 
independently by means of dispersion theory. This formalism allows one to deduce the 
universal pion-pion final-state interactions from the accurately known tttt phase shifts and, 
in the scalar sector, a coupled-channel treatment with the kaon-antikaon system. In this 
work an analysis of the decays J/'ip7r~^TT~ and is presented. We 

find very good agreement with the data up to 1.05 GeV in the Trvr invariant mass, with a 
number of parameters reduced significantly compared to a phenomenological analysis. In 
addition, the phases of the amplitudes are correct by construction, a crucial feature for 
many CP violation measurements in heavy-meson decays. 
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1 Introdnction 

H-meson decays can be exploited for Standard Model tests and beyond, in particular to 
determine the Cabibbo-Kobayashi-Maskawa (CKM) couplings and to study CP violation. 
For a theoretical description of many of these decays, it is mandatory to understand the 
strong final-state interactions in terms of amplitude analysis techniques [1], with tight con¬ 
trol over the magnitudes and phase motions of the various partial waves involved. For 
example, the decays B —)• /o(980)K5 and B —)• (j){1020)Ks are explored for an experi¬ 
mental determination of the CP asymmetry sin 2/3 [2-5], /3 being one of the angles of the 
unitarity triangle, which requires precise knowledge of the strange and non-strange scalar 
form factors that we discuss in this article. We focus on the decays B^ — J/'ip7r~^7r~ and 
Bg —>■ J/'07r+7r“, measured by the LHCb collaboration [6, 7]. The tree-level process of the 
weak decay into J/if) and a qq pair is depicted in Fig. 1 (exemplarily for the B^ decay). 
These analyses complement former related studies of B^ and B^ decays by the BaBar [8], 
Belle [9], CDF [10], and DO [11] Collaborations as well as older LHCb results [12, 13]. 
Universality of final-state interactions dictates that the hadronization into pions and the 
rescattering effects in the 7r'’'7r“ system for S- and P-waves are closely related to the scalar 
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Figure 1. The — >■ J/'07r+7r“ diagram to leading order via W~ exchange. The hadronization 
into pions (S'-wave dominated) proceeds through the pion strange scalar form factor T® (s). In the 
case of the —>■ J/i/itt+tt” decay, with sod, the pions are generated out of a non-strange scalar 
source, i.e. T® (s) is replaced by the pion non-strange scalar form factor r))(s) for S'-wave and by 
the vector form factor for P-wave pions. 

and vector pion form factors, respectively. We describe these form factors using dispersion 
theory, using Omnes (or MuskhelishvilPOmnes) representations. In doing so we exploit 
the fact that LHCb found no obvious structures in the invariant mass distribution, 

suggesting that left-hand-cut contributions in the vr+vr" system due to the crossed-channel 
J/^7r+ interaction are small and can be neglected. 

The advantage of the dispersive framework is that all constraints imposed by analyt- 
icity (i.e., causality) and unitarity (probability conservation) are fulfilled by construction. 
Further, it is a model-independent approach, so we do not have to specify any contributing 
resonances or conceivable non-resonant backgrounds. For the vector form factor a single¬ 
channel (elastic) treatment works very well below 1 GeV. In the scalar sector the strong 
coupling of two S'-wave pions to KK near 1 GeV due to the /o(980) resonance, causing a 
sharp onset of the KK inelasticity, necessitates a coupled-channel treatment. Therefore a 
two-channel Muskhelishvili-Omnes problem is solved. This two-channel approach breaks 
down at energies where inelasticities caused by 47r states become important, we are thus 
not able to cover the complete phase space, but restrict ourselves to the low-energy range 
y/s < 1.05 GeV. 

In Ref. [6] the decay is described by six resonances in the tt+tt" channel, /o(500), 
/9(770), a;(782), /?(1450), p(1700), and /2(1270), which are modeled by Breit-Wigner func¬ 
tions. This parametrization of especially the /o(500) meson is somewhat precarious, as 
the broad bump structure of this scalar resonance is not well described by a Breit-Wigner 
shape. As demonstrated for the first time in the context of B decays in Ref. [14], it 
should be replaced by the corresponding scalar form factor. In the present work this idea 
is extended and rigorously applied using form factors derived from dispersion theory. In 
particular, there is no need to parametrize any resonance, since the input required to 
describe the final-state interactions is taken from known phase shifts, and therefore the 
/o(500) appears naturally in the non-strange scalar form factor. The B^ decay, described 
in the experimental analysis by hve resonances, /o(980), /o(1500), /o(1790), /2(1270), and 
/2(1525) (Solution I) or with an additional non-resonant contribution (Solution II), dom¬ 
inantly occurs in an 5-wave state [7], while the P-wave is shown to be negligible. Given 
the almost pure ss source the pions are generated from, this decay shows great promise to 
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provide insight into the strange scalar form factor. 

The idea of such a “scalar-source model”, where an 5-wave pion pair is generated out 
of a quark-antiquark pair and the final-state interactions are described by the scalar form 
factor, is also used in Ref. [15] for the description of the Bg and decays into the scalar 
resonances /o(980) and /o(500), respectively. It was employed earlier e.g. in analyses of 
the decay of the J/ij: into a vector meson (a; or cf)) and a pair of pseudoscalars (vrvr or 
KK) [16, 17]. In these references the strong-interaction part is described by a chiral 
unitary theory including coupled channels, which yields a dynamical generation of the 
scalar mesons. In contrast to the present study, the very precise information available on 
pion-pion [18-21] and pion-kaon [22] phase shifts is not strictly implemented there. Related 
studies using the chiral unitary approach are performed in Ref. [23], where the J/?/;-vector- 
meson final state is analyzed, and in Ref. [24], which includes resonances beyond 1 GeV. In 
contrast to models of dynamical resonance generation, the scalar resonances are considered 
as qq or tetraquark states in Ref. [25]. Other theoretical approaches employ light-cone QCD 
sum rules to describe the form factors [26] . Progress on the short-distance level is made in 
Ref. [27], where the factorization formulae (which we treat in a naive way) are improved 
in a perturbative QCD framework. 

This manuscript is organized as follows. In Sec. 2, we review the construction of the 
transversity amplitudes and partial waves, after sketching the kinematics. We provide ex¬ 
plicit expressions that relate the theoretical quantities to the angular moments determined 
in experiment. Section 3 is focused on the Omnes formalism. The fits to the LHCb data, 
using the —>• J/'i/’vr^vr" angular moment distributions, are discussed in Sec. 4, where 

we use several configurations with and without D-wave corrections to study the impact 
of certain corrections to our fits. We also predict the 5-wave amplitude for the related 
Bg JK~ decay. The paper ends with a summary and an outlook in Sec. 5. Some 
technical details are relegated to the appendices. 

2 Kinematics, decay rate, and angular moments 

In this section we derive the decay rate and angular moments for the B^ —?• 

decay mode in terms of partial-wave amplitudes up to D-waves, employing the transversity 

formalism of Ref. [28] . The formalism works analogously for the Bg decay. 

2.1 Kinematics 

The kinematics of the decay B'^^^{pb) «7/'!/'(p^)vr+(pi)7r“(p2) {J/'4’ A*"*"!*”) can be 

described by four variables: 

• the invariant dimeson mass squared, s = {pi +P 2 )‘^, 
and three helicity angles, see Fig. 2, 

• Gj/^, the angle between the in the J/V’ rest frame and the J/ij: in the B^jg 

rest frame (S^); 
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Figure 2. Definition of the kinematical variables for - 8 °/^ —>■ J. 


• 07 r, the angle between the vr^ in the vr^vr center-of-mass frame and the dipion 
line-of-flight in 

• (f), the angle between the dipion and the dimuon planes, where the latter originate 
from the decay of the J/V'. 

The three-momenta of either of the two pions in the dipion center-of-mass system (p^r) and 
of the J/V’ in the rest frame (p^) are given by 


2 ^/s 2 ’ 2 m B ms 

with the Kallen function A(a, b, c) = + b'^ + — 2ab — 2ac — 26c. 

We define the two remaining Mandelstam variables as 


= —> ( 2 - 1 ) 


( 2 . 2 ) 


t = {PB- Pif and u = {pB- P2f, 
where the difference of these two determines the scattering angle 0 .^-, 

t - u = - 2 p^{pi - P 2 ) = - 2 aT,X cos 0 ^. 

Further, we introduce two additional vectors as combinations of the above four-momenta, 

P^=Pi+P2^ Q^=Pi-P2- ( 2 - 4 ) 


(2.3) 


2.2 Matrix element 

To calculate the matrix element we make use of the effective Hamiltonian that governs the 
6 —>■ ccd transition [29] , 

= ^{VM [Ci{p)0^{p) + C2{p)02{p)] + ...}, (2.5) 

where the Ci are Wilson coefficients and the Oi local current-current operators 

Oi = Ck 7 f,(l - 75)6/ di7^(l - j 5 )ck = Ck 7 f,(l - 75)cfc d/7^(l - 75)6^, 

O 2 = Ck 7 ^(l - 75)6fc dz7^(l - 75 )q = c^ 7^,(1 - 75)01 d/7^(l - 75 ) 6 ^, ( 2 . 6 ) 
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with k, I being color indices. In the second step the quark operators are regrouped by means 
of a Fierz rearrangement. The ellipses in T-LeS denote operators beyond tree-level, including 
penguin topologies. Vcb and Vcd are the CKM matrix elements for c ^ b and c ^ d (where 
Vcd is to be replaced by Vcs for the decay), and Gp = 1.166365 x 10“® GeV“^ is the 
Fermi constant. 

Under the assumption that the final-state interaction between the J/V’ and the pions 
is negligible (no obvious structures are found in the J/V’tt channel experimentally [6, 7], 
and a close-to-zero J/ipn scattering length aj/^^ = —0.01(1) fm results from lattice cal¬ 
culations [30]) a factorization approach appears to be justified. Note that on the quark 
level this naive factorization ansatz may be spoiled [31, 32], for instance due to (large) 
penguin contributions that we have neglected in Eq. (2.5) [33, 34]. However, a more com¬ 
plicated structure of the source term does not conflict with our approach; any factorization 
limitations due to color structures do not concern the hadronic hnal-state interaction, for 
which the short-distance factorizations are sufficient but not mandatory. All we use is the 
fact that the B decays provide clean qq sources of much shorter range than that of the 
final-state interaction. In our approach, any deviations from clean point sources would be 
parametrized by derivatives of the source term. An excellent fit to the data even without 
those correction terms is a proof that with respect to the final-state interactions the sources 
can be regarded as point-like. 

We express the matrix elements of the four-quark operators by two independent hadronic 
currents, valid if the cc system produced by the hadronization of the virtual W~ is well 
separated from the spectator quark system. For the decay of the B^ meson considered here 
the matrix element is, in analogy to the B^ expression given in Ref. [35], written as 

7r+(pi)7r"(p2)|d7M(l - 75 )^ 1 ^°/^(ps)) X {J/il;{p^,e)\c'y>^c\0), 

-.- ^ - ' '---' 

-Mr = (2.7) 


Mf, = 


with o®® = Ci{p)+C 2 {n)/Nc-\-. ■ ■, the ellipses denoting combinations of Wilson coefficients 
due to penguin diagrams, we have not taken into account explicitly. 

The scale (/r) dependence of the Wilson coefficients is cancelled by the scale dependence 
of the hadronic matrix elements, cf. Sec. 3; p is chosen to be of order 0{mB), such that 
heavier particles, in particular the W, are integrated out. 

The current that creates the J/'ijj from the vacuum is related to the decay constant f^. 
The matrix element containing the pions is given by the three transversity form factors Wq, 
W||, and J-±, corresponding to the orthogonal basis of momentum vectors [28] 


a -a 

Pp Pi±) = 




A 


-iP^p)aPl3Q'f, Qn\) =Q'"- 


{P ■ Pi^){Q ■ pu s{Q-p^) 


A2 


X2 P^' 


( 2 . 8 ) 


We define e 


fii/pa 


such that €0123 = —= -|-1. 
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The partial-wave expansions of the transversity form factors read^ 

= ^v/27TTT-r(s)P^(cos0^) 

e 


= (s) -I- a/S cos (s) + ^ (3 cos^ - 1 ) (s) -I-... , 




V2ITT 


VW+^) 




)P^'(cos 6»^) = yi" (s) -F V5 cos (s)^ -h • • •, 

(2.9) 


where the ellipses denote waves larger than D-waves. In Appendix A the relation to the 
helicity form factors is briefly sketched, which have a well-known partial-wave expansion. 


2.3 Decay rate and angular moments 

When comparing the angular moments to the experimental data we have to deal with flavor- 
averaged expressions due to the B^-B^ mixing and take into account the CP-conjugated 
amplitudes (the B^ decay mode) as well. Since the interfering term between the amplitudes 
is negligibly small [6], the decay rate can be written as the sum of the decay rates for the 
direct B^ and the mixed CP-conjugated B^ mode, 


d^r (P^ ^ J/'^TT+vr ) d^r (direct) d^T (P° ^ ) 

d d cos Ott dy^ d cos dy^ d cos 0.^ 


( 2 . 10 ) 


Note that this neglect is less justihed when applying the formulae to the B^ decay rate. In 
the analysis of Ref. [7] an interference term is added to Eq. (2.10). However, in Sec. 4.2 
we find that it is sufficient to take into account 5-waves. In that case the interference 
term does not affect the fit procedure and merely generates a tiny shift of the resulting fit 
parameter (the normalization Cq). 

In this section we provide expressions for one particular mode. The CP-related ampli¬ 
tude can be deduced straightforwardly by multiplying the transversity partial-wave ampli¬ 
tudes with CP eigenvalues as outlined in detail below (cf. the discussion around Eq. (2.15)). 

The differential decay rate is given by 


d^T 




dy/i d cos 0 T, 


4(47r)3m 


3 yyi ^ 
B 


(s) + \/3 cos (s) + ^ (3 COS^ 0jr - l) (s) 


+ sin^ ^ (s) + + VE cos (s) ^ 


( 2 . 11 ) 


^ Though we expect the D- and higher waves to be small and therefore describe only S- and P-waves in 
the Omnes formalism, we present the formulae including the D-wsve contribution, as we will study their 
impact at a later stage. 
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see Appendix A for details. By weighting this decay rate by spherical harmonic functions 
yj°(cos 6 * 7 r), we define the angular moments 

. ll°(cosg^)dcosg^. (2.12) 

y_i dysdcos 0 -^ 

With the orthogonality property 

/■i X. . 

/ yj°(cOs6*7r)K°(cOS07r)dcOS07r = —, (2-13) 

J-i 27r 


we obtain 



V^{Y 2 ^) 




2(47r)3m| 


B 


+ a. 

'i 2 


(P) 


+ 




(p) 


+ 




T 

•J f 


{D) 


{S) 


+ 


T 

•J r 


(p) 


+ 


F 


(D) 




(47r)3m| 


+ 




B 


Re ( 


Re ( 


+ Re ( F^Pf!:^'^* 




(47r)^m| 


B 


Re ( fPfP‘ 


+ 


F 

■J r 


(D) 


^Re[FPFP' 


P f 

2^5 V 


f: 


{p) 


+ 


F 


(P) 


+ 


Mr 

14 V 


F 


(D) 


+ 


+ 


V5 


I 

■J r 


(P) 


F 


(D) 


+ 


V5 


T 

•J r 


(D) 


(2.14) 


where (Iq^) corresponds to the event distribution, (Rj^) describes the interference between 
S- and P-wave as well as P- and D-wave amplitudes, and (Y^) contains P-wave, P-wave, 
and S~D-wave interference contributions. 

The corresponding expressions for the CP-conjugated modes are related to the above 
equations by certain sign changes due to the CP eigenvalues rjcp = ±1 in the definitions 
of the transversity partial-wave amplitudes, as already mentioned in the beginning of this 
section. We declare the amplitudes Fr^'^ to describe the decay, then the corresponding 
decay amplitudes are given by 



TlCpX\ 


{(.) 


(2.15) 


with rjcp = -|-1 for the r = 0, || P-waves and the r =T D-wave, and rjcp = —1 otherwise. 
Consequently the angular moments (Rq^) and {Y 2 ) unchanged under CP conjugation, 
while the conjugated moment (Rj^) has opposite sign, such that when considering flavor- 
averaged quantities and summing over the B^ and B^ contributions, (Rj®) vanishes. In the 
following we thus consider (Rg^) and {Y 2 ) only. 


3 Omnes formalism 

We describe the S- and P-wave amplitudes using dispersion theory. This approach allows 
us to treat the pion-pion rescattering effects in a model-independent way, based on the 
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fundamental principles of unitarity and analyticity: the partial waves are analytic functions 
in the whole s-plane except for a branch-cut structure dictated by unitarity. In the following 
we deal with the functions // (s) (referring to isospin / and angular momentum £) that 
possess a right-hand cut starting at the pion-pion threshold sthr = 4M^ and are analytic 
elsewhere, i.e. we do not consider any left-hand-cut or pole structure related to crossing 
symmetry. This is justihed from the observation that there are practically no structures 
observed for the crossed J/'i/’vr"’' channel in the region of interest [6]. 

Considering two-pion intermediate states only, Watson’s theorem holds, i.e. the phase 
of the partial wave is given by the elastic pion-pion phase shift [36], and the discontinuity 
across the cut can be written as 

disc//(s) = f/(s-hie} - fi(s-ie) = 2ia^f/(s) [4(s)]* = //(s)e"*'^«%indf. (3.1) 

A solution of this unitarity relation can be constructed analytically, setting (compare 
Ref. [37]) 

f^{s) = P{s)ni{s), (3.2) 

where P{s) is a polynomial not fixed by unitarity, and the Omnes function Q^(s) is entirely 
determined by the phase shift Sj(s) [38], 

Oj(.) = .xp|ir dA, (3.3) 

s'(s' - s-t€) j 

with 

n^(0) = 1 and ^jis) / 0 Vs. (3.4) 

The P-wave amplitudes can be well described in the elastic approximation up to energies 
of roughly 1 GeV.^ The simplest possible application is the pion vector form factor (s), 

(0|j^m(0)k’^(7’i)7r"(P2)) = {P 2 - PiTF^(s), (3.5) 

which obeys a representation like (3.2) with a linear polynomial Ppv{s) = 1 -|- as, a ~ 
0.1 GeV“^ [39] up to yfs « 1 GeV, with the exception of a small energy region around 
the uj resonance that couples to the two-pion channel via isospin-violating interactions. In 
this context it is important to note that the electromagnetic current jem, introduced in 
Eq. (3.5), can be decomposed as 

iem = \ {ui^u - dy'^d) i (uy^u -h dy'^d). (3.6) 

Thus it contains with the first term an isovector and with the second term an isoscalar 
component. The latter couples directly to the oj, whose decay into tt+tt” is suppressed 
by isospin, but enhanced by a small energy denominator (i.e., the small width of the u), 


^ In the following we will suppress the isospin indices as Bose symmetry demands the 5-waves to be 
isoscalar, while the P-waves are restricted to 7 = 1. 
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hence leading to a clearly observable effect in the pion form factor [40-42], Theoretically, 
this effect is correctly taken into account by the replacement [43-45] 

PFv{s)n\{s) Ppv{s)nlis)(^i + ^2 -s)- 

Note that in case of the u the use of a Breit-Wigner parametrization is appropriate since 
the CO pole is located far above the relevant decay thresholds and since = 8.5 MeV is 
very small. A fit of the form factor parametrization introduced in Eq. (3.7) to the KLOE 
data [41] yields Kem ~ 1-8 x 10“^. This fixes the strength of the so-called p-co mixing 
amplitude phenomenologically. The isospin-violating coupling Kem is of the usual size, 
however, near the co peak its smallness is balanced by the factor ~ 90 from the co 

propagator, giving rise to an isospin-violating correction as large as 15% on the amplitude 
level, corresponding to 30% in observables due to interference with the leading term. Note 
also that the p-co mixing amplitude has been pointed out to significantly enhance certain 
CP-violating asymmetries in hadronic B-meson decays [46]. 

The effect of the co on the —>■ J/tt~ decay can be related straightforwardly to that 
on the pion vector form factor. To see this observe that the source term for the tttt system 
is dd at tree level, see Eig. 1, such that the isospin decomposition of the corresponding 
vector current reads 

d-f^^d = - ^ ^ + dy'^d). (3.8) 

Comparison to Eq. (3.6) shows that the relative strength of the isoscalar component differs 
from the electromagnetic current by a factor of —3, such that we will fix the p-co mixing 
contribution in analogy to Eq. (3.7), but with the replacement Kem —>■«: = —3Kem ~ 
—5.4 X 10“^. Notice that this is in contrast with the experimental analysis [6], where the 
CO contribution is fitted with free coupling constants. 

The (elastic) single-channel treatment, introduced in the beginning of this section, can¬ 
not be used in the 5-wave case: there are strong inelastic effects in the region around 1 GeV 
due to the opening of the KK channel, coinciding with the /o(980) resonance, which affects 
the phase of the scalar pion form factors (see e.g. the discussion in Ref. [47]). Thus the 
Omnes problem has to be generalized, with the Watson theorem fulfilled in the elastic re¬ 
gion and inelastic effects included above the KK threshold. This leads to the two-channel 
Muskhelishvili-Omnes equations that intertwine the pion and kaon form factors, defined 
as 


(7r+(pi)7r (p 2 ) \qq\ O) = R'^r^(s), 

{K+{p,)K-{p 2 )\qq\ 0> = R'?r^(s), (3.9) 

where the quark flavors may be either qq = {uu+dd)/2 for the light quarks, with the super¬ 
script q = n denoting the corresponding scalar form factor, or qq = ss for strange quarks 
(with superscript q = s). Eurthermore, = M^/{mu + rrid), R® = (2M|. — M^)/{2ms). 
Note that the form factors r^^(s) are invariant under the QCD renormalization group, 
while the hadronic matrix elements are not due to the scale dependence inherent in the 
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factors This in turn allows for the cancellation of the scale dependence in the Wilson 
coefficients introduced in the effective Hamiltonian of Sec. 2. 

Appealing to the tree-level diagram of Fig. 1, we expect the non-strange scalar form fac¬ 
tors to contribute dominantly in the decay, while the strange ones should feature mainly 
in the corresponding decay of the Bg. As discussed in detail below, these expectations are 
confirmed by the data analysis. 

The Muskhelishvili-Omnes formalism is briefly reviewed in Appendix B. It requires 
three input functions: in addition to the vr-Tr phase shift already necessary in the elastic 
case, modulus and phase of the tttt ^ KK 5-wave amplitude also need to be known. 
Our main solution is based on the Roy equation analysis by the Bern group [20, 21] for 
the TTTT phase shift, the modulus of the tttt —>■ KK 5-wave as obtained from the solution 
of Roy-Steiner equations for ttK scattering performed in Orsay [22], and its phase from 
partial-wave analyses [48, 49]. Alternatively, we employ the T-matrix constructed by Dai 
and Pennington (DP) in Ref. [50]: here, a coupled-channel R-matrix parametrization is 
fitted to rnr data [51-55], and the Madrid-Krakow Roy-equation analysis [19] is used as 
input; furthermore, the KK threshold region is improved by fitting also to Dalitz plot 
analyses of Df [56] and D+ — K^K~tt^ [57] by the BaBar Collaboration. 

In addition, the channel coupling manifests itself through the fact that even in the sim¬ 
plest case, corresponding to the polynomial of Eq. (3.2) reducing to a constant, the scalar 
form factors depend on two such constants, corresponding to the form factor normaliza¬ 
tions for both pion and kaon. In contrast to the single-channel case, here the shape of the 
resulting form factors depends on the relative size of these two normalization constants; on 
the other hand, once this relative strength is fixed, it relates the final states yttf and KK to 
each other unambiguously. We will make use of this additional predictiveness in Sec. 4.3. 

In order to apply this formalism to the transversity partial waves we have to construct 
partial waves /r^^(s) that are free of kinematical singularities, i.e. represented by functions 
whose only non-analytic behavior is related to unitarity. In Appendix A the hadronic 
matrix element is introduced (using the basis of the momenta p^, P^, and Q^, Eq. (2.4)) 
in terms of the form factors Ai and V, Eq. (A.l), and related to the transversity basis, 
Eq. (A. 5). Given that the form factors Ai and V are regular, Eq. (A. 5) implies that there 
are additional factors of X, and ^/s introduced into the transversity form factors, which 
give rise to artihcial branch cuts in the unphysical region. To avoid those, we write the 
partial waves as 

T-f )(s) = X/f )(s), J-f )(s) = )(s), 

.fP(s) = PPis) = ^Xf[^\s), (3.10) 

where the are treated in the Omnes formalism, i.e. 

fPis) = P^^'^\s)r:{s) + P^^’^\s)rUs), /f)(5) = PPis)Q\{s). (3.11) 

Eor the 5-wave, we a priori allow for contributions of both non-strange (n) and strange 
(s) scalar form factors. The coefficients of the polynomials Pt^\s} are to be determined 
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from a fit to the efficiency-corrected and background-subtracted LHCb data, in particular 
to the angular moments {Yq) and 

Basically we assume the various polynomials to be well approximated by constants. 
However, to study the impact of a linear correction at a later stage, we also consider linear 
polynomials = 6”(1 -|- 6 g”s) and Pr^'^ = ar{l -|- a'^s) for the non-strange S-wave and 

the P-wave amplitudes, respectively. The strange S-wave contribution is expected to be 
very small (in the LHCb analysis of —)• J/'ip7r~^TT~ the /o(980) meson is not seen), but 
tested in the fits. On the contrary, the Bg J/'ip7r~^TT~ distribution is dominated by the 
/o(980) resonance, described by a constant polynomial times Omnes function, = Cq, 

while there is no structure in the /o(500) region reported by LHCb. Thus in that case the 
non-strange S-wave amplitude is assumed to be negligible, to be confirmed in the fits. 

Although the first P-wave resonance seen is the /2(1270), it may affect also the region 
below ^/s « 1 GeV due to its hnite width, T= 185.1^2 4 MeV [58]. Therefore we also 
test its influence on the ht. The P-waves could be treated in the same dispersive way 
as S- and P-waves, but this would increase the number of free parameters in our fits to 
the LHCb data. As the effect of P-wave corrections is rather small, we avoid introducing 
additional fit parameters and take over the amplitudes (with hxed couplings) used in the 
LHCb analysis, where the /2(1270) resonance is modeled by a Breit-Wigner shape. 

Since the data are given in arbitrary units, we collect all prefactors in normalizations 
that we subsume into the fit parameters (and into the transversity coefficients ap that we 
extract from the LHCb fit results). Writing (Y^) in terms of Omnes functions for S- and 
P-waves, supplemented by the P-wave resonance contribution, yields 



V4^(h'2°) 


Xa^V~s!^X^\b^{l + b',^s)T^is) + cgr^(s)|' 

-I- al |hJ(s)|^ ([ao(l + “qs)]^ + 's[a||(l + aj|s)]^ -|- sX^[a_L(l -|- a'_i_s)]^^ 

T=0,±,|| ' 


Xa^V~s |2Re + 6[r5)r)[(s) + 4K{s)\ (s) 

2 

-I- |hJ(s)|^ (^ 2 [ao(l -I- Aqs)]^ - 's[a||(l + Oys)]^ - sX^ [a±(l -b 


-b 


V5 

V5 


4 ( 0 ) 




(3.12) 


For details concerning the definition of the Breit-Wigner amplitudes r = 0, ||,T, 

see Ref. [ 6 ]. 


4 Fits to the LHCb data 
4.1 —)• J 

We fit the angular moments (Lg') and (T 2 °), Eq. (3.12), simultaneously. Taking up the 
discussion of Sec. 3, our basic ht, FIT I, includes three ht parameters (to be compared 


- 11 - 










to 14 free parameters in the Breit-Wigner parametrization used in the LHCb analysis, 

(P) 

see below): the normalization factors for the S-wme (6g) and for two P-waves /q and 
/|l (ao,a||). (We find that including the r =± P-wave amplitude practically does not 

change the i.e. a_L is a redundant parameter.) In the basic fit only S- and P-waves 
are considered. Beyond that, we study the relevance of certain corrections: in FIT II we 
use again the same three parameters as in FIT I, but in addition we include the P-wave 
contributions, fixed to their strengths as determined by LHCb. To further improve FIT II, 
supplemental linear terms (6g,aQ,aj|—cf. Eq. (3.12)) are allowed in FIT III. Performing 
FIT III we find that two of the slope parameters, the linear non-strange 5-wave term (6 q) 
and the r =|| P-wave slope (aj|), yield no significant improvement of the fits; their values 
are compatible with zero within uncertainties. We thus fix them to zero, and in FIT III only 
the four parameters 6g, oq, ay, and Og are varied. Furthermore, the effect of an inclusion 
of a strange 5-wave component is tested. Its strength is found to be compatible with zero, 
justifying its omission. 

Note that the scalar pion form factors depend on the normalizations of both the pion and 
kaon form factors. While the normalizations in the case of the pion form factor are known 
quite precisely, there are considerable uncertainties for the kaon form factor normalizations, 
having an impact on the shapes of both pion form factors, see Appendix B. The non-strange 
kaon normalization F^(0) is limited to the range (0.4... 0.6). In our fits we fix the value 
to r^(0) = 0.5, which is compatible with the current algebra result. The effect from a 
variation of F^ (0) in the allowed interval shows up only in the second decimal place of the 
X^/ndf. 

The fitted coefficients and the resulting x^/ndf, referring to Eq. (3.12), are listed in 
Table 1. The large uncertainties can be traced back to the correlations between the fit 
parameters, especially present in FIT III. For a comparison to the LHCb fit, we insert 
their fit results (best model) into our definition of the x^- more specific terms this 
means that we do not compare to the x^ published in Ref. [6], for which the full energy 
range up to ^/s = 2.1 GeV is fitted with 34 parameters and the data of all angular moments 
(Y^) for i = 0,..., 5 are included, but we calculate the x^ in the region we use in our fits, 
i.e. including data up to ^/s = 1.02 GeV and the angular moments (kg') and (1^2°) only. 
We obtain XLHCb/'^hf = 2.08. In this limited energy range the Breit-Wigner description, 
including the /o(500), p{770) and t<;(782), requires 14 fit constants, while we have three 
(FIT I, H) or four (FIT HI) free parameters and find x^/ndf = 2.0 (FIT I), x^/ndf = 1.5 
(FIT H) and x^/ndf =1.3 (FIT HI). The calculated angular moments for the three fit 
models in comparison to the data are shown in Fig. 3. 

Probably the most striking feature of our solution is the pronounced effect of the uj that 
leads to the higher peak in Fig. 3. As mentioned above, this isospin-violating contribution 
is fixed completely from an analysis of the pion vector form factor, however, its appearance 
here is utterly different, since the coupling strength is multiplied by a factor of —3. This 
not only enhances the impact of the co on the amplitude level to about 50%, but also 
implies that the change in phase of the signal is visible a lot more clearly: while in case 
of the vector form factor the a; amplitude leads to an enhancement on the p-peak and 
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Figure 3. (Yq) (left) simultaneously fitted with {Y 2 ) (right), using 3 parameters without _D-wave 
contribution (FIT I, red, solid), and improving step by step by adding a Breit-Wigner-parametrized 
D-wave contribution (FIT II, blue, dashed) and by allowing for 4 free parameters, also supplemented 
by the D-wave contribution (FIT III, green, dotted). 



xVndf 


ko 

kii 

FIT I 

1.97 


46.5l®:° 

51.81?,- 

FIT II 

1.54 


-0-6.6 

49.51?,- 

FIT III 

1.32 


o'y 'y-|-20.3 
'^*•*-21.3 

0-4lo:7 


Table 1. Resulting fit parameters and x^/ndf for the various fit configurations FIT I-III for the 

—>■ J/'ijjTT^Tr~ decay. 

some depletion on the right wing, forming a moderate distortion of the line shape, here 
we obtain a depletion on the p-peak accompanied by an enhancement on the right wing. 
While the current data do not show the oj peak clearly, a small shape variation due to the 
p-u interference is better seen in Ref. [33], where a finer binning is used. The p-uj mixing 
strength obtained from a fit in that reference is consistent with the strength we obtain in a 
parameter-free manner. Nonetheless, improved experimental data are called for, since an 
experimental confirmation of the uj effect on —)• J would allow one to establish 
that the decay indeed provides a rather clean dd source. 

A key feature of the formalism employed here is its correct description of the S-wave. 
Figure 4 shows the comparison of the 5-wave amplitude strength of the LHCb Breit- 
Wigner parametrization with the ones obtained in FIT I-III, as well as the comparison of 
the corresponding phases. In the elastic region, the phase of the non-strange scalar form 
factor dp" = arg(r”) coincides with the tttt phase shift dg that we use as input for the 
Omnes matrix, in accordance with Watson’s theorem. Right above the KK threshold, dp'* 
drops quickly, which causes the dip in the region of the /o(980), visible in the modulus 
of the amplitudes as well as the non-Breit-Wigner bump structure in the /o(500) region. 
We hnd that the phase due to a Breit-Wigner parametrization largely differs from the 
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Figure 4. Comparison of the S'-wave amplitude strength and phase obtained in the LHCb and 
in our hts, respectively. In the left panel the S'-wave part of the decay rate for the three fit 
configurations FIT I-III is depicted together with the LHCb outcome. The right panel shows the 
phases of the non-strange scalar form factor dp" (equal to the tttt S-wave phase shift dg below the 
KK threshold) compared to the S-wave phase d/g extracted from the LHCb analysis. 

dispersive solution, indicating that parametrizations of such kind are not well suited for 
studies of CP violation in heavy-meson decays. 

Note that in the analysis of Ref. [33] the /o(500) is modeled not by a Breit-Wigner 
function, but by the theoretically better motivated parametrization of Ref. [59]. In this 
work, higher resonances are included by multiplying S-matrix elements. While this proce¬ 
dure preserves unitarity, it produces terms at odds with any microscopic description of the 
coupled irn-KK system. As such also this approach introduces uncontrolled theoretical 
uncertainties into the analysis. The only stringently model-independent way to include 
hadronic final-state interactions is via dispersion theory. 

4.2 —)• J/'i/>7r+7r~ 

The Bg J/'^TT+TT” distribution in the region up to roughly 1 GeV is clearly dominated 
by the /o(980). We therefore describe the data with the strange S'-wave component only, 
using a constant subtraction polynomial (cq). The only non-zero contribution to the fit thus 
comes from (To°)- Fitting the data up to ^/s = 1.05 (1.02) GeV yields x^/ndf = 2.2 (1.8) 
and Cq = 16.8±0.4 (16.8±0.4). In analogy to the decay we also perform the fit including 
the D-wme parametrization of the LHGb analysis. This yields an additional non-zero 
contribution to (T^^) due to the S-D-wave interference, which is fitted simultaneously with 
(Tq^). Further, the influence of a linear subtraction polynomial for the strange S-wave is 
tested. However, none of these corrections exhibits a considerable improvement. 

In the LHCb analysis the full energy range, ^/s < 2.1 GeV, is fitted with 22 (24) 
parameters for Solution I (H). Confining to the region we examine in our fit and considering 
the /o(980) resonance only, the number of fit parameters reduces to four (six), and we 
calculate ~ (0-82), when using our definition of the x^- 


- 14 











Figure 5. Left panel: (Yq^) fitted using the strange S'-wave with constant subtraction polynomial 
for two different phase inputs (red, solid: B+0 input [20-22], green, dotted: DP input [50], based on 
the Madrid-Krakow analysis [19]). Right panel: comparison of the phase of the strange scalar pion 
form factor for the B+0 (blue, dashed) and DP (red, solid) input, respectively, with the S'-wave 
phase extracted from the LHCb analysis (Solution I and II, shown with error bands). 


The strange scalar form factor, or the /o(980) peak in the dispersive formalism, depends 
crucially on the Trvr KK 5-wave transition amplitude, which is not as accurately known 
as elastic TTTr scattering (and even contains subtleties as non-negligible isospin breaking 
effects due to the different thresholds of charged and neutral kaons, see e.g. Ref. [60]). As 
there are no error bands available for the Omnes matrix (or the various input quantities), 
to estimate the theoretical uncertainty we use and compare the fits resulting from the two 
different coupled-channel T-matrices described in Sec. 3. A minimization of the using 
the modified Omnes solution based on Ref. [50] yields x^/ndf = 3.4 (2.4) and Cq = 18.3 ± 
0.5 (18.2 ± 0.5).^ The resulting (Tg^) curves for both fits, using the phase input from the 
Bern [20, 21] and Orsay [22] groups (B+0), as well the one of Ref. [50] (DP), are presented 
in Fig. 5. Furthermore we show the phase shifts and the phases of the strange form factor 
for both phase inputs and compare to the LHCb phase due to Solution II (with /o(980) 
and a non-resonant 5-wave contribution) as well as Solution I (/o(980) parametrization 
only). While the latter phase has a negative slope for s < 1 GeV, which does not agree 
with the known phase shift, the phase extracted in Solution II is remarkably close to both 
the Bern and Madrid phase motions. 

4.3 ^ J/ijjK+K- S -wave prediction 

Having obtained the —>■ fit parameters, we can make a prediction for the 

JK~ 5-wave amplitudes, using the relation between the Trvr and the KK final 
states provided by the coupled-channel formalism, cf. Appendix B.^ 

® A similar procedure for the 13° decay has a rather small effect since the S-wave is not dominant in that 
case, and the difference of the P-wave phase of Refs. [19-21] is quite small (the S- or P-wme phase 
modification yields, in the most perceptible cases, a 4% correction of the y^). 

■* In the case of the 13° ^ J/'tpK'^K~ decay [61], the prediction of the S'-wave does not work in such a 
direct way due to the 7 = 1 S'-wave contribution (with a prominent ao(980) resonance) in addition to fo 
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In particular an understanding of the 5-wave background to the prominent (/>(1020) is of 
interest. In the LHCb analysis [62], the /o(980) as well as a non-resonant 5-wave content 
is reported within a mass window of ±12 MeV around the (/>(1020), which contribute an 
5-wave fraction of (1.1 ± 0.1 ]*^q^)%— consistent with former measurements from LHCb, 
CDF, and ATLAS [63-65], as well as theoretical estimates [66]. We calculate the 5-wave 
fraction in the same mass interval ±12 MeV around the (^(1020) mass adopting the LHCb 
Breit-Wigner parametrization for the i;/)(1020), but using the predicted 5-wave for the 
J/il^K^K~ final state. Naively, this 5-wave can be obtained by replacing the pion scalar 
form factor and all pion masses and momenta by the respective kaon quantities and taking 
the resulting fit parameters from the pion case. However, the fit result depends on the 
normalization of the —)• J/'^vr+vr" distribution. Hence, taking over the pion fit results 
for such a prediction requires a proper normalization of both decay channels relative to 
each other. To achieve this, we use the absolute branching fractions [58] 

B (5° ^ J/i^K+R-) = (7.9 ± 0.7) x 10"^, 

B ^ J/V^TT+Tr") = (2.12 ± 0.19) x 10"^ 


and define normalization constants 

{^R} N{B0 ^ J/2lj{Tr+7r-,K+K-})' 

where 

V(H° ^ J/iP{Tr+TT-,K^K-}) = j (Vo” (5° ^ J/V'{7r+7r", AT+A"})) d^s (4.2) 
is the total number of events.^ 

The 5-wave contribution to the 0(1020) peak region is given by 


T^s/4> 


C-Vllv (fl» ^ JI,I,K+K-) ) aR 


(4.3) 


where we can approximate the (normalized) angular moment in the region of interest by 
the 5-wave and the 0(1020) contribution, 


\\/~s—m^\<12 MeV 


(4° ^ J/0A+A-) 

« XaKV~s[x^M^ |cgr|,(V^)|' ±Wx |«M5r^(s) 

' T 


(4.4) 


Using the B±0 input, we obtain 'R,s/(t> — 1-1%, in agreement with the LHCb result. How¬ 
ever, there is a notable uncertainty due to the estimated ambiguity in the phase input in 
the region of the /o(980) resonance discussed in Sec. 4.2. Using the DP phase instead of 
the B±0 phase input yields a fraction of 1.95%. 


resonances in the 1 = 0 S-wecve. 

® For the J/ipK'^ K~ decay [62] no data for the efficiency-corrected angular moments are available. 

We therefore extract the strength of the 0(1020) Breit-Wigner amplitude from the published expected 
signal yield Nexp and use N(B° J/0N+N-) = B{B° -)■ J/ipK+R-)/N,^p. 
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5 Summary and outlook 

In this article, we have described the strong-interaction part of the J/'k~ and 

—>■ J/^7r+7r“ decays by means of dispersively constructed scalar and vector pion form 
factors. This formalism respects all constraints from analyticity and unitarity. The non- 
strange and strange scalar form factors are calculated from a two-channel Muskhelishvili- 
Omnes formalism that requires the pion-pion elastic S'-wave phase shift as well as modulus 
and phase of the corresponding tt-tt —?■ KK amplitude as input. For the vector form 
factor, an elastic Omnes representation based solely on the pion-pion P-wave phase shift 
is sufficient, supplemented by an enhanced isospin-breaking contribution of p-ui mixing, 
which can be fixed from data on e'’'e“ ^ tt+tt”. 

For energies ^/s < 1.02 GeV, a minimal description of all S- and P-waves (constructed 
in a form free of kinematical singularities) as the corresponding form factors, multiplied 
by real constants, has been shown to be sufficient. Allowing for subtraction polynomials 
with linear s-dependence leads to a slightly improved fit quality solely in the case of one 
P-wave component, with a slope still compatible with zero within uncertainties. In partic¬ 
ular considering the 5-wave slope as a free fit parameter (as opposed to fixing it to zero) 
only yields a minimal improvement of the In accordance with expectations from the 
underlying tree-level decay mechanism, below the onset of P-wave contributions that be¬ 
come important with the /2(1270), only the non-strange scalar and the vector form factors 
feature in the decay, while the strange scalar form factor determines the P® 5-wave. 

The overall fit quality in the energy range considered is at least as good as in the 
phenomenological fits by the LHCb collaboration [6, 7], where Breit-Wigner resonances 
and non-resonant background terms were used. However, since the dispersive analysis 
allows one to use input from other sources, our analysis calls for a much smaller number 
of parameters to be determined from the data. In addition, a comparison of the P^ 5- 
wave obtained from the dispersive analysis with the one deduced from the LHCb analysis 
shows drastic differences in both modulus and phase: it is well-known that the /o(500) does 
not have a Breit-Wigner shape, and therefore such parametrizations should be avoided— 
especially when it comes to studies of CP violation that need a reliable treatment of the 
phases induced by the hadronic final-state interactions [14]. The LHCb analysis of the 
P° 5-wave uses a Flatte parametrization of the /o(980), solely (corresponding to their 
Solution I) or combined with a non-resonant background (Solution H). Only Solution H 
yields a phase that is close to the phase of the strange scalar form factor, and approximately 
compatible with Watson’s final-state interaction theorem in the elastic region. 

Finally we have made a prediction for the P^ —)• J/xl:K~^K~ 5-wave, which is related 
to the corresponding 7r'*“7r“ final state through channel coupling. Only the results of the 
fit to the 7r'’'7r“ final state are required to predict an 5-wave fraction below the i;i)(1020) 
resonance of about 1.1%, in agreement with the findings by the LHCb collaboration. We 
have not attempted a corresponding prediction for the P^ —)• J/'ipK~^K~ 5-wave, since 
this has an isovector component (corresponding e.g. to the ao(980) resonance). This would 
have to be described by a coupled-channel treatment of the vrr/ and KK 5-waves [67]. 

To extend our description of the form factors to higher energies, eventually covering most 
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of the energy range accessible in —)• J/tt~ , inelastic channels with corresponding 

higher resonances have to be taken into account. Here, a formalism recently developed for 
the vector form factor [45] that correctly implements the analytic structure and unitarity, 
reduces to the Omnes representation in the elastic regime, but maps smoothly onto an 
isobar-model picture at higher energies should be extended to the scalar sector. Even 
an extraction of the scalar form factors from these high-precision LHCb data sets seems 
feasible, and should be pursued in the future. 
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A Form factors and partial-wave expansion 


In the standard basis of momenta p^, P^, and Q^, Eq. (2.4), the matrix element describing 
the hadronic part of the decay is given by four dimensionless form factors, three axial 
{Ai) and one vector (V), 

= (vr+(pi)7r"(p2)| = P^Ai + Q^A2 + ip^)^iA2, + 

Jl = (A.l) 

In Sec. 2.2 we use a different (orthogonal) basis of momentum vectors, and 

see Eq. (2.8), corresponding to the orthonormal basis of polarization vectors of the 
J/^ meson [28], 


(*) = #. '"(O) =-^To)' 


1 


M, 




cjjr sin 0 ^, 


(QqI) T (^-2) 


This allows us to describe the matrix element in terms of the transversity form factors, 
Eq. (2.7), or similarly (with regard to an easily performable partial-wave expansion) in 
terms of helicity form factors, defined via the contraction of with the polarization 

vector, 


P, = (7r+7r-lJ^"/*lB,^/,)eT(A). (A.3) 


The relations between the transversity and helicity form factors can be read off to be 


Pt = At, Bo = Ao, P± = (A|| ± A^)^ sin^^e^'-^, 


(A.4) 
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as well as those to the set {Ai, V} 


J-± = —^/sXV, J^ll = ■\/sA2, Xq = 


X 


2A, + 

X 


P-Q . (T^X cos I 

•'t — ~rj —- 4.1 + 


M,, 


2XM^ 


{{P-p^){P • Q) - sMl) A 2 + M^As. 


(A.5) 


The unphysical time component Pt does not contribute. We expand the remaining three 
form factors Po,± in partial waves. The latter relation is of particular interest when defining 
partial waves that are free of kinematical singularities and zeros, see Sec. 3. 

The partial-wave expansion of the helicity amplitudes reads 


-Hxis) = ^ VW+Tnf\s)dio{9^)e^^^, (A.6) 

i 


where the are the small Wigner-d functions. Using 


doo(^7r) — P(.{pOs6A): diQ(d,r) — ~d^ig(07r) — — 


sin 


^/W+^ 


Plicos 9^), 


(A.7) 


we see that the zero-component T-Lo{s) is expanded in terms of Legendre polynomials 
Pi{cos 9^) and thus contains all S-, P-, and D-wme contributions, while the 'H±{s) partial- 
wave expansions, proceeding in derivatives of the Legendre polynomials ^^(cosd^), start 
with the P-wave amplitudes, i.e. 


Po{s) = Pq^\s) + V3cos9,,Ti^ff\s) + ^ (3cos^ - l) Pq^\s) , 

^^±(■ 5 ) = T-^^sind^ -I-\/5cosd7r7^llf^('S)) -b ..., (A.8) 


where the ellipses denote P-waves and larger. Equivalently, due to Eq. (A.4) and using 
(s) = {^) d=-7'j^^(s)), we arrive at the partial-wave expansion of the transver- 

sity form factors given in Eq. (2.9) in the main text. 

In order to calculate the differential decay rate we sum over the squared helicity ampli¬ 
tudes, 

\M\" = (|PoP + \n+\^ + \n.\^) {q = {d, s}) (A.9) 

and integrate over the invariant three-particle phase space, which is given by 


d$(3) 


AcTtf 

4(47r)2m^ 


ds d cos djr d^. 


(A.IO) 


Neglecting waves larger than P-waves and integrating over cj) we arrive at Eq. (2.11). 


B Coupled-channel Omnes formalism 

We briefly discuss the coupled-channel derivation of the scalar pion and kaon form factors 
(/ = 0,.^ = 0). The two-channel unitarity relation reads 

discr(s) = 2iT^*{s)^{s)V{s), (B.l) 
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where the two-dimensional vector r(s) contains the pion and kaon scalar isoscalar form 
factors and Tq{s) and S(s) are two-dimensional matrices, 



2icTj,(s) 


7yO(s)e2d^o(«)-'5o(«)) 

2i<j^(s) 


\ 

/ 


(B.2) 


and E(s) = diag (cr7r(s)0(s — 4M^), crK(s)0(s — 4M^)), with (Tj(s) = (l — and 

0(.) denoting the Heaviside function. There are three input functions entering the T- 
matrix, the tttt S'-wave isoscalar phase shift 5o(s) and the tttt —)• KK S-wave amplitude 
5 (q(s) = |5 (Q(s)|exp(z'0Q(s)) with modulus and phase. The modulus |5o(s)| is related to the 
inelasticity parameter r/o(s) by 


ho(s) = ^1 - 4o-^(s)cJi^(s)|ff[}(s)|20(s - 4M|). (B.3) 

Writing the two-dimensional dispersion integral over the discontinuity (B.l) leads to a 
system of coupled Muskhelishvili-Omnes equations, 


r 



T0*(s')S(s')r(s') ^^, 
s' — s — ie 


(B.4) 


A solution can be constructed introducing a two-dimensional 0mnes matrix, which is 
connected to the form factors by means of a multiplication with a vector containing the 
normalizations r 7 r( 0 ) and ri^(O) [68], 


/ T^{s) ^ / r^(0) \ 

V J [ n2lis) ^22(8) ) ^r^(o) ) ’ 

where r,r,ic(s) represents both strange and non-strange form factors, T® ^(s) and T” ^(s), 
which differ merely in their respective normalizations. Thus the problem reduces to finding 
a matrix Q(s) that fulfills 


1 poo J^0*( „r\y( 

imn{s) = T^*{s)j:{s)n{s), n{s) = - o 01 01 q(o) = 1, (b. 6 ) 

vr s' - s - le 


which has to be solved numerically [68-71]. To ensure an adequate asymptotic behavior, 
we exploit the correlation between the high-energy behavior of the 0mnes solution and the 
sum of the eigen phase shifts ^ h^is) [69], 




rriTT 


TT for / = 1, £ = 1 J s^oo 1 

n^[s) —-, 

27r for / = 0, = 0 ^ 


(B.7) 


where m is the number of channels that are treated in the formalism. 

According to the Feynman-Hellmann theorem, the form factors for zero momentum are 
related to the corresponding Goldstone boson masses, which at next-to-leading order in the 
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Figure 6. Modulus of the scalar pion non-strange (left panel) and strange (right panel) form factors, 
depicted for three different normalizations inside the allowed range, illustrated by the uncertainty 
band. 

chiral expansion in terms of quark masses depend on certain low-energy constants. These 
are determined in lattice simulations with Nf = 2 + 1 dynamical flavors at a running scale 
// = 770MeV [72], limiting the form factor normalizations to the ranges® 

r;((o) = 0.984 ± 0.006, r^(o) = (- 0 . 001 ... 0 . 006 ) « 0 , 

r^(0) = (0.4...0.6), rf^(O) = (0.95... 1.15). (B.8) 

Figure 6 shows the results obtained for the modulus of the pion form factor (see also 
Ref. [74]). The sensitivity due to the uncertainty in the kaon form factor normalization 
is illustrated by the uncertainty bands. The strange form factor exhibits a peak around 
1 GeV, which is produced by the /o(980) resonance. On the contrary in the pion non- 
strange form factor the a meson appears as a broad bump (notice the non-Breit-Wigner 
shape) around 500 MeV. 
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